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Abstract 

The trace anomaly of the energy-momentum tensor is usually quoted in the 
form which is proportional to the beta function of the theory. However, there 
are in general many definitions of gauge couplings depending on renormalization 
schemes, and hence many beta functions. In particular, N = 1 super symmetric 
pure Yang-Mills has the holomorphic gauge coupling whose beta function is one- 
loop exact, and the canonical gauge coupling whose beta function is given by the 
Novikov-Shifman-Vainshtein-Zakharov beta function. In this paper, we study which 
beta function should appear in the trace anomaly in J\f = 1 pure Yang-Mills. We 
calculate the trace anomaly by employing the M = 4 regularization of M = 1 pure 
Yang-Mills. It is shown that the trace anomaly is given by one- loop exact form if the 
composite operator appearing in the trace anomaly is renormalized in a preferred 
way. This result gives the simplest resolution to the anomaly puzzle in J\f = 1 pure 
Yang-Mills. The most important point is to examine in which scheme the quantum 
action principle is valid, which is crucial in the derivation of the trace anomaly. 



1 Introduction 



In supersymmetric field theories, any operator is contained in some sup er mult iplet. This 
is also the case for the energy-momentum tensor T^ v '. It was discovered long time ago [T] 
that the energy-momentum tensor is contained in a supermult iplet called the supercurrent 
. The components of the supercurrent is given as 

J* = f R + OS" + h.c. + 29a u 9T^ + ■■■, (1) 

where is an R-symmetry current and is the current of supersymmetry. 

In the past, it was believed [2] that the conservation equation of the energy-momentum 
tensor is extended to the superspace equation given by 

D & J a a « D a X (2) 

where X is some chiral superfield, or else 

T>°Ja& ~ ~\D 2 D a J, (3) 

where J is some real vector superfield, and we use the symbol ~ for equations which are 
valid if equations of motion are used (i.e., the equations are valid up to contact terms 
when they are inserted into correlation functions). Eq. ([3]) is possible only if a conserved 
R-symmetry current exists, so we focus on Eq. in the following. 

Eq. (j2J) gives the conservation equations of T^ v and S 1 *, but it also gives the additional 
constraint given by 

l^ + id^^Fjc, (4) 

where Fx is the F-component of X. The fact that the right hand side of Eq. (jlj) is given 
by a single F-term causes a problem called the anomaly puzzle [3]. The axial anomaly 
contributing to c^j'r is often said to be one-loop exact due to the Adler-Bardeen theo- 
rem [4] . This one- loop exactness seems to indicate that the chiral superfield X is one- loop 
exact. On the other hand, the trace anomaly usually receives higher order corrections. 
Then, higher order corrections to X seem to be necessary. These two statements look 
inconsistent with each other. This is the anomaly puzzle. Many works have been done to 
solve the anomaly puzzle, including Refs. [51 EJ [7J |8] and references therein. 



2 



The situation have been changed by the discovery of a new supercurrent multiplet [9] . 
(See also Ref. [10] for an early work.) It has been found that the supercurrent equation 
can be relaxed to the form 

D«J a a w D a X - -D 2 D a J. (5) 



Then, Eq. (jJJ) is now modified to 



\ T H - MFx)- 1 - 

d,f R « Im(F x ), (7) 



-T; « Re(F x )--L> J; (6) 



where Dj is the D-component of J. These equations suggest that X can be one-loop exact 
to maintain the one-loop exactness of the axial anomaly, and higher order corrections to 
the trace anomaly is accounted for by the new term Dj. Based on this observation, 
a resolution to the anomaly puzzle has been proposed [IT] . (See also Ref. [12] for a 
related work.) More explicitly, let us consider a theory with matter chiral fields $ r which 
transform in the representation r of the gauge group. The gauge field strength chiral field 
is denoted as W a , and the superpotential is denoted as W. Then, X and J are given as 

4 

.XII - \ <I> 

d$ r 32vr 



X 

3 



- E ^ - 3tAd ~f tr w*w a 



(8) 



J = -4j2lMe- 2V ®r, (9) 

r 

where tjM, and t r are the Dynkin indices for the adjoint representation and the represen- 
tation r, respectively, 7 r is the matter anomalous dimension, and the summation is over 
all matter fields. Several consistency checks have been done for this proposal in Ref. [11] . 

However, there still remains a puzzle. The problem already exists in M = 1 supersym- 
metric pure Yang-Mills, and we restrict our attention to this case. In this theory, there 
is no candidate for the operator J, and hence the situation regarding the anomaly puzzle 
of this theory have not been changed. The trace anomaly is often quoted in the form 



where (3 is the beta function of the gauge coupling, and F^ u is the gauge field strength. 
In supersymmetric theories, there is a known "exact" beta function [131 E], called the 
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Novikov-Shifman-Vainshtein-Zakharov (NSVZ) beta function. In the case of the M = 1 
pure Yang-Mills, it is given as 

dg 2 1 3t Ad g* 

/W - ^ - -^2 1 _ tAdg 2 /8n 2- C") 

If this beta function is inserted into Eq. ({TO]) , the trace anomaly is not one- loop exact. 
Then it is inconsistent with the Adler-Bardeen theorem and the supercurrent equation. 

However, the beta function is obviously renormalization scheme dependent. For exam- 
ple, the beta function in the so-called DR (or DR) scheme is not given by the NSVZ beta 
function beyond the two- loop level [13] . We can even define a gauge coupling, called the 
holomorphic gauge coupling, whose beta function is one- loop exact [7j. Furthermore, the 
composite operator F )JLV F* IV is also renormalization scheme dependent. Eq. ( ITOj) has been 
proven [151 HI EE] by using the minimal subtraction scheme in dimensional regularization 
for QCD-like theories, but this scheme cannot preserve manifest supersymmetry. Because 
the anomaly puzzle originates from the superspace equation, we have to use a scheme in 
which supersymmetry is preserved in a manifest way. 

In this paper, we perform calculations of the trace anomaly by employing the M = 4 
regularization of the Af = 1 pure Yang-Mills [H [18], which is suitable for our purpose. 
(Regularization by dimensional reduction (DRED) [191 120] can also maintain manifest 
supersymmetry^ But this scheme is not suitable for the study of the anomaly puzzle, 
as we explain in section [5j) We will find that the trace anomaly is indeed given by the 
one-loop exact form, if the operator F^F^ is renormalized in a way which is preferred 
by supersymmetry and the topological nature of F^F^, where F^ v = \e^ up<T F^. 

In section |2l we give an indirect derivation of the trace anomaly in general field theories, 
and clarify the point at which the scheme dependence enters. In section [3] the M = 4 
regularization is introduced. In section HI the trace anomaly formula is applied to the 
M = 4 regularized version of the Af = 1 theory. There we discuss how the scheme 
dependence becomes important. Finally, conclusions are given in section [5j We also 



discuss the relation between our work and the previous works [5j El [Tl [8] . Appendix A 
contains a review of basic properties of the energy-momentum tensor. 



: We will not care about mathematical inconsistency, violation of supersymmetry, and problems with 
75 in DRED. See Ref. [21] and references therein for discussions on these things. 
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2 Trace anomaly formula 



In this section, we give an indirect derivation of the trace anomaly (up to a total derivative 
term) in a general field theory which has a Lagrangian description. We will clarify the 
assumption which is made in the derivation of the trace anomaly. 

Let us consider a theory described by a set of fields (fi, where / is a label specifying 
the fields. We assume that the theory is renormalized in some renormalization scheme. 
We consider a correlation function, 

(n^ite)) = / [DtplexpiiSfaKVUtpi&i), (12) 

where A a are parameters (e.g., couplings and masses) of the theory, and S is the action. 
We assume that the field </?/ has mass dimension Dj and the parameter A a has mass 
dimension d a . By a simple dimensional analysis, we obtain 

= ("E + D,) + 4 + £ 4A. A) , (13) 

where fj, is the renormalization scale. On the other hand, the renormalization group (RG) 
equation tells us that 

where 7/ is the anomalous dimension of (fi, and j5 a is the beta function defined by 

d 

(3a = V-q-K- (15) 

We assume that there are no operator mixings among (pi for simplicity. 

Next, by using the Ward-Takahashi identity of the energy-momentum tensor given by 
Eq. ( 1A.7j) of Appendix A we obtain 



Ij ^T;(y)n^o)=E(^^+ d ^) (n^(^))- (is) 

where d is the space-time dimension, and Cj are arbitrary constants representing the am- 
biguity in the definition of the energy-momentum tensor, (see Appendix A for details). 
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We have assumed that the integral of the total derivative / d d xd fJ ,(x u T> lu ) vanishesE Com- 
bining Eqs. ffT3l) . (TH|) and ffl~6]) . we obtain 

viM)h ( 17 ) 



where we have defined Aj = Dj + 7/. We can use Eq. flA.9[) to rewrite the first term of 
the right-hand-side as 

(y: {dc h - a 7 j) =-«'((/ ^e m - a /) ^/(y)) n^o*)) > 

where iVj is an operator which vanishes by equations of motion (see Appendix A). 



Now we make a crucial assumption. Naively, by using the path integral expression for 
the correlation function, the derivative with respect to A a is given as 

, n<P*(xi)) = (i[d<v^ll<pi i {x i )), (19) 



d\ a W 7 \ J d\ 

where C is the Lagrangian. This is known as the quantum action principle 
Assuming that this is really the case, we obtain 



d d yT»(y)Hvu 



(/ dd y (E - dCi) N M + EK^ - /y ^7) • (20) 

Therefore, we finally obtain the operator equation for the trace anomaly, 

r)C. 

TjH = £(d«A« - 0a)^- + E ( A ' - dCi) Nj + d^j v , (21) 

a OA a j 

where d^j v is a total derivative term which we do not try to determine. As noted above, 
Nj vanishes by equations of motion. Therefore, as expected, vanishes up to the total 
derivative term in the fixed point of RG flow, d a X a — (3 a = 0. 

In the above derivation, the use of the quantum action principle (Tl9l) is crucially im- 
portant. The validity of the quantum action principle depends on what regularization or 



2 We should subtract the "cosmological constant term" in T^ u so that the vacuum expectation value 
is given by (T^) = 0. 
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renormalization scheme is used. For example, it is known that Eq. (1191) holds when min- 
imal subtraction in dimensional regularization is used to renormalize all the parameters 
and operators of the theory [23J . Then, Eq. ( 12T|) reproduces the result of Refs. [151 [H3 E] 
for the trace anomaly in QCD-like theories up to a total derivative term. 

If Eq. ffl9l) is naively applied in the Af = 1 pure- Yang- Mills, we would obtain 

<} ((■■■)) = i (( [d*xd 2 d]w a W a + h.c) (■■■)) (22) 



dg- 2 NV " \\J 4 

where (• • ■) represents a chain of operators. If this is the case, the trace anomaly is 
proportional to (j3/g 4 )W a W a . In the following sections, we examine in what scheme 
Eq. ( 12 2 p is really valid in the Af = 1 pure Yang-Mills. 

3 Af = 4 regularization of Af = 1 pure Yang-Mills 

We employ the Af = 4 regularization of the Af = 1 pure Yang-Mills which is discussed in 
Refs. [H [IE] . In this regularization, we take the Af = 4 theory with Af = 4 supersymmetry 
softly broken to Af = 1 by the mass terms of the adjoint chiral fields. The Lagrangian is 



given by 



c = [ d 2 e-^w a w a + h.c. + f d 2 ed 2 eY,zMe- 2V ^i 



1=1 



+ J d 2 6 {y/2f ABC ^M + E + h.c, (23) 

where A, B and C are gauge indices, /abc is the structure constant of the gauge group, 
and gauge indices are omitted in other terms. The parameters g\ and Zi (i = 1, 2, 3) are 
given in terms of the gauge coupling g and the theta angle 9 as 



Z X Z 2 Z Z = Re(^ 2 ), (25) 



9? = 9- 2 -^ (24) 



where the second equation comes from the constraint of Af = 4 supersymmetry. Following 
Ref. |18j, we have normalized the fields so that the coefficient of the term fABC&i&2 



3 We follow the notation and convention of Wess and Bagger [21] , except that a matter kinetic term 
is given by $^e _2y $, and a gauge field strength chiral field is given by W a — ^D 2 (e 2V D a e~ 2V ). 
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does not depend on the gauge coupling g. In this normalization, we can maintain the 
holomorphy [25J when the holomorphic coupling g^ and the mass terms Mj (i = 1,2,3) 
are extended to background chiral superfields. 

In the energy region much below the scale Mj, this theory becomes the same as the 
M = 1 pure Yang-Mills, up to terms which are power suppressed by Mj. However, the 
existence of the chiral multiplets <3>j makes the theory UV-finite thanks to the known 
fmiteness of the N = 4 theory^ Therefore, this theory can be regarded as the J\f = 1 
pure Yang-Mills regulated by the adjoint chiral fields with the cutoff scale Mj. 

Before closing this section, some technical remarks are in order. Just to make our 
arguments more concrete, we further regulate the theory by e.g., regularization by dimen- 
sional reduction (DRED). One reason that we introduce the further regularization is that 
we would like the regularization to be consistent with the extension of parameters into 
background superfields so that we can utilize the power of holomorphy. Strictly speak- 
ing, this extension maintains only M = 1 supersymmetry, and a possibility exists that 
the M = 4 regularization might not be compatible with the extension of couplings into 
background superfields. However, we believe that there is no such incompatibility. For 
example, let us suppose that the couplings g^ 2 and Z\ (and also Mj) are allowed to have 
non-vanishing constant F-terms (and D-terms for Zj) with the constraint ( )25|) imposed as 
a superfield equation. Then, the violation of M = 4 supersymmetry is only soft, i.e., it is 
broken only by parameters with positive mass dimensions. Therefore all the divergences 
in dimensionless quantities are still cancelled, and we need no counterterms for the lowest 
components of g^ 2 and Z*. Other soft breaking parameters are contained in the same 
multiplet as g^ 2 or Zj, so it is expected that no counterterms are required also for these 
soft breaking parameters. This idea may be elegantly realized in the scheme of analytic 
continuation of parameters into superspace discussed in Ref. [26]. For example, RG equa- 
tions have been given for soft breaking terms by analytically continuing the RG equations 
for g^ 2 and Zj into superspace, and hence all the soft breaking terms are manifestly RG 
invariant if that is the case for the lowest components. This suggests that all the UV 
divergences are cancelled in the present theory even if the parameters are analytically 

4 It is important that the J\f = 4 supersymmetry is only softly broken by the mass terms of the adjoint 
chiral fields, and hence the mass terms do not affect the UV divergences of dimensionless quantities. 
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continued as in Ref. (26] fl 

An important check of the cancellation of divergences in our regularization framework 
can be found in Ref. [26] . When g~^ 2 and Zj have nontrivial soft terms, a counterterm for 
the mass of the so-called epsilon scalar in DRED is required. At the one-loop level, this 
counterterm, when analytically continued into superspace, is given as 



- / d 2 9d 2 9- 1 ' 



16^ 



t Ad \og(Re(g h 2 )) - U lo § Z i 



(I^rO, (26) 



where fi is the Lorentz index of the "compactified" 2e dimensions in DRED, tAd and tj are 
Dinkin indices for the adjoint representation and the representation of the matter field % 
respectively, and Y ^ is defined as 

r„ = \afD^ v D a e- 2V ). (27) 

One can check that is gauge covariant, and hence the above term is allowed by gauge 
invariance [27]. If Eq. ( 126]) were present, this would give a contribution to the gauge 
kinetic term by using D 2 {Y^Y0 = 4eW a W a . In our theory, we have tj = tAd (i — lj 2, 3), 
so Eq. ( 1261) vanishes thanks to the condition ( l25i) . For this cancellation to happen, it 
is important that the adjoint chiral fields are normalized in the way which preserves 
the holomorphy of the superpotential [18], leading to Eq. f )25|) . We expect that this 
cancellation continues to happen at higher order level. In the following, we assume that 
the M = 4 regularization is compatible with the extension of the parameters to M = 1 
superfields (at least if they are space-time constants, i.e., d^g^ 2 = etc.), and the usual 
argument based on holomorphy is justified in our regularization framework^ 

Another reason we regularize the J\f = 4 theory by e.g. DRED is that we need to treat 
composite operators. Even in M = 4 theories, composite operators require regularization 
and renormalization in general. However, fortunately, this problem is also not essential 
for our purpose of computing the trace anomaly. In our analysis, the only composite 



5 Here we have introduced the framework of Ref. [26] not because we want to have soft breaking 
parameters in our theory, but because we want to make the statements about holomorphy more explicit. 
For our purposes, it is important that we are working in a well-defined (i.e., explicitly regularized) setting 
where the argument of holomorphy can (at least in principle) be made rigorous. 

6 Much easier way to ensure finiteness may be to use AT = 2 supersymmetry, by taking Z\ — Re(<?^ ) 
and Z% = Z% = 1. In this case, we can extend <?^ 2 to a background TV = 2 vector multiplct [28,. Then 
the absence of radiative corrections is established. 
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operators which we treat are the energy-momentum tensor and the operators appearing 
in the trace anomaly. The energy-momentum tensor is finite (at least if improvement [29] 
is done appropriately). On the other hand, we need the operator d£/dg~ 2 for the trace 
anomaly. For simplicity, let us neglect the mass terms for the adjoint chiral fields. Then 
this operator is indeed finite in M = 4 theories. The most easy way to see this is to note 
that the coupling g is exactly marginal, i.e., the coupling does not run in the RG flow for 
all the possible values of it. Then the operator dC/dg^ 2 should have scaling dimension 4 
(at least up to possible total derivative terms which we do not care), which coincides with 
its mass dimension. Therefore there is no anomalous dimension for this operator. This 
fact indicates that no renormalization is required for this operator. 

4 The trace anomaly in J\f = 1 pure Yang-Mills 

Now we discuss the calculations of the beta functions and the trace anomaly in the M = 1 
pure Yang-Mills. Let us first discuss the beta functions. In the previous section, we have 
argued that the manifest holomorphy is expected to be maintained in the regularization 
process. Then, the low energy physics depends only on the combination 

8tt 2 3 

— log A/i. (28) 

9h i=\ 

This statement can be confirmed directly at the one-loop level (by decoupling the adjoint 
chiral fields, perhaps in the manifestly supersymmetric framework of Ref. [26]). Higher 
order corrections are absent due to the holomorphy and the dependence of g^ 2 on the 
theta angleE 

We introduce a single cutoff scale M (which is taken to be a real positive number) and 
will give the masses Mi in terms of M. We have infinitely many choices for doing this [H] , 




This argument for the one-loop exactness becomes invalid when we try to apply it for SQCD reg- 
ularized by finite Af = 2 theories. In the low-energy theory, we have to regard matter wave-function 
renormalization as an independent parameter in addition to the holomorphic coupling, as is evident in 
the analyses of e.g. Ref. |26) . We have two (or more) parameters in the low-energy theory, but there 
is only one gauge coupling g in the high-energy theory. Then it is not possible to choose g so as to fix 
both the holomorphic gauge coupling and the wave-function renormalization of the low-energy theory 
when the cutoff scale is changed. This difficulty can be avoided by rescaling the matter fields so that the 
wave function renormalization of the low energy theory becomes unity. However, this process violates 
the holomorphy, and hence the one-loop exactness is lost. 
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and this choice determines the renormalization prescription for the bare coupling g. We 
will study the RG equation for g which makes the combination (1281) invariant when the 
cutoff scale M is changed. Although we discuss the RG equation for the bare quantity, 
it is also easy to define renormalized couplings in our regularization framework. We will 
comment on this point later. 

Two particular choices for Mj are often discussed. One choice is given by Mj = M, 
which maintains the holomorphy about the cutoff M. In this case, by requiring that the 
combination ( )28|) is invariant under the change of M, we obtain 

M mf - < 29 > 

Thus we obtain the one-loop exact beta function in this scheme. The other choice is that 
we take Mj = ZjM, so that the cutoff M corresponds to the tree level masses of $j when 
these fields are canonically normalized. Then, the real part of Eq. (12811 becomes 



% ~ E tAd log(^M) = ^--t Ad log g- 2 - 3t M log(M) (30) 

where Eq. ( I25p has been used. This choice leads to the NSVZ beta function 

as d 2 1 3t Ad g* 

M m 9 = (31) 

Now we would like to discuss the trace anomaly. Before doing that, a remark on the op- 
erator W a W a is necessary. The trace anomaly is given in terms of the composite operator 
F lxv F tJ,v , which is contained in W a W a . This superfield also contains the operator F^F^. 
Due to its topological property, the operator F^ U F^ V should not receive multiplicative 
renormalization, and hence the operator W a W a also should have a preferred overall scale. 
(However, see section [5] for discussion on trouble in DRED if the M = 4 regularization is 
not used.) This does not mean that the operator F iW F iW requires no renormalization at 
all [201 EI]- Although the overall scale of F^F^ should be determined by its topological 
nature, radiative corrections of the form F llv F^ v — > F PjV F^ v + d^j^ is not restricted by 
topological arguments, where is a gauge invariant operator. In the language of differen- 
tial forms, we can add an exact 4-form d(*j) to the closed 4-form F A F without affecting 
it as an element of the de Rham cohomology group. For example, in the case of QED, the 
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Figure 1: An example of the divergence of the operator F^F^ . Wavy lines are gauge 
fields, and solid lines are fermions. The black filled circle is the operator F^F^ . 

one-loop diagram in Figure [TJ gives a divergence and hence the operator F I1U F^ V needs a 
counterterm given by 

[V1 = F, U F^ + ( -JL-t + c finitc J d,jl (32) 

where we have used dimensional regularizationjj] j£ = §^[7^,75]^ is the axial current, 
and [F^F^] is the renormalized operator. The term Cfi n i te is a finite counterterm whose 
value depends on renormalization prescriptions, and hence the operator [F^F^] has an 
ambiguity of the form c Rnite d^j^. 

One might think that the distinction between F I1V F^ U and is meaningless (in 
massless limit of fermions) because the axial anomaly equates them. If that were the 
case, F^F^ would get multiplicative renormalizations. However, the operator d^j^ is 
equivalent to F flu F fJ ' u only if we use equations of motion, i.e., d^ji? ~ F^ U F^ V /8n 2 , and they 
differ by the presence of contact terms when they are inserted into correlation functions. 
The contact terms are essential in the discussion of the topological properties of F^F^. 
Meaning of the Adler-Bardeen theorem also requires careful considerations. See Ref. [TT] 
for discussions on the importance of this observation in solving the anomaly puzzle. 

In supersymmetric theories, the radiative corrections and counterterms discussed above 
are extended to superspace as 

[W a W a ] = W a W a + D 2 J } (33) 

where J is a real vector superfield. Therefore this type of corrections should be allowed 

in the definition of the composite operator W a W a . 

8 Here we have avoided DRED because DRED is extremely subtle to perform this type of calculations, 
because of the relation tr(7 5 7 M 7 I/ 7 p 7 cr ) = (see e.g. Ref [H]). We hope that this 75 problem does not 
cause trouble as long as we are working in superspace without looking components of superfields. 
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Let us calculate the trace anomaly using the formula (l21j) . In the previous section, 
we have regularized the theory further by DRED for concreteness. It is known that the 
quantum action principle (|T9|) holds for bare quantities in DRED [21]. Thus we can use 
Eq. (1191) in our regularization. The important point is that we have to include regulator 
contributions which come from terms containing $j. 

If we use the scheme Mj = M which have led to the one-loop exact beta function, we 
obtain § 

dC r „„1 



dg- 



2 



J d 2 6- (w a W Q + D 2 K) + h.c. (34) 

As explained in the previous section, the operator d£/dg~ 2 is finite due to M = 4 su- 
persymmetry. However, the operator W a W a itself receives divergent corrections even in 
J\f = 4 theories from diagrams such as Figure [TJ Thus the term D 2 K can be regarded 
as a regulator term which makes the operator W a W a finite. We define the renormalized 
operator as|^ 

[W a W a ] R = W a W a + D 2 K. (36) 

As explained above, the addition of D 2 K does not violate the topological nature of F^F^, 
and hence it should be allowed. Therefore, the trace anomaly is given as 



3tAd 



J d 2 9[W a W a } R + h.c. +drf$ (37) 



327T 2 

where m means that we have used the equations of motion, Nj pa 0. Therefore, we have 
determined that the coefficient of the operator [W 70 !^]^ in the trace anomaly is one- loop 
exact. 



9 We also have to include contributions from gauge fixing terms and ghosts. However, we expect that 
these contributions are BRST exact, as in the case of QCD |T5] [TSJ [T7] . We neglect them in this paper. 

10 The operator K, and hence [H^ a Wa]jj, may look to be depending on the precise choice of Zi, but it 
is not so. If we use different wave function renormalizations Z[ satisfying Eq. (|25[) and obtain K 1 , then 
the difference K' — K is a current of SU(3) C SU(4)# symmetry of the Af = 4 theory. In particular, it 
is UV-finite. Then the argument of the heavy field decoupling theorem (see e.g. Ref. [32]) may tell us 
that K' — K is power suppressed by the heavy field masses M,. The UV-finiteness is essential, because 
otherwise momentum integrals such as J d 4 k(k 2 + M 2 )~ 2 spoil the suppression. 
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Next, let us perform the calculation in the scheme Mj = ZiM, which have led to the 
NSVZ beta function. Then, we obtain 

^ = /^(l™ fi+ ig^M i$ f) + h.c. (38) 

In this case, we have additional contributions coming from the regulator mass terms. 
Therefore, the naive use of Eq. (122j) is invalid if the regulator contribution is not taken 
into account seriously. In the low energy limit (i.e., in the energy region much below M), 
the operator M^f becomes 

^ -£rz[W a W a ] R (2 = 1,2,3). (39) 

At the one-loop level, this equation can be checked directly. The holomorphy restricts 
higher order corrections to the form W a W a — > W a W a + D 2 J. We have assumed that these 
corrections are such that the operator W a W a is replaced by [W a W a ]ji. This is because 
the operator is UV-finite since it is protected by M = 4 superconformal symmetry. 
Using Eq. (1391) . we obtain 

8C f j2 nlfi tAd^dfogZ, 



d 2 e] i- t -^-)[w a w a } R + h. c ., (40) 



4 V' 8?r 2 j 

where we have used Eq. (1251) . Combining this result with the NSVZ beta function ( 13 ip . 
we again get the same answer given by Eq. (l37j) . Therefore, the one- loop exact form of 
the trace anomaly is also confirmed in this scheme. 

Finally, let us comment on renormalization in the M = 4 regularization. We have 
treated bare couplings up to now, but it is very easy to define a renormalized coupling 
since we know that the low energy physics depends only on the combination (128]) . For 
example, we can define a renormalized coupling g R as 

where \x is the renormalization scale. Now we can let the cutoff Mj be infinity at the end of 
the calculations with g R fixed. The derivation of the trace anomaly is completely parallel 
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to the one given above. Although it is not so evident that the operator [W a Wo]fl remains 
finite in the limit Mj — > oo, we expect that this is the case since the energy- momentum 
tensor (and hence its trace) is UV-finite. 

5 Conclusions and comparison with the literature 

In this paper we have investigated the trace anomaly in M = 1 pure Yang-Mills to solve 
the anomaly puzzle in this theory. The trace anomaly is often quoted as — (fS/Ag^F^F^, 
but the validity of this formula depends on whether the quantum action principle given 
by Eq. (|22l) (or more generally Eq. ( Fl9l) ) is really true or not. To settle this issue, we have 
explicitly regularized the theory by using the J\f = 4 regularization. Then the regulator 
contribution to Eq. ( 1221) can be explicitly seen, and this contribution depends on what 
scheme we use. In the scheme in which we obtain the one-loop exact beta function, it 
is almost correct to use Eq. f )22|) naively with the regulator contribution neglected. The 
regulator contribution only plays a role of making the operator W a W a renormalized, 
[^"^Jk = W a W a + D 2 K, without spoiling the topological nature of F^F^ contained 
in W a W a . However, if the scheme is used in which we obtain the NSVZ beta function, 
the regulator contribution is not at all negligible, as seen in Eq. ( 1401) . If this contribution 
is correctly taken into account, we obtain the same answer for the trace anomaly in both 
of the schemes. The result for the trace anomaly has the one-loop exact form, and hence 
is consistent with the supercurrent equation and the Adler-Bardeen theorem. 

Although the M = 4 regularization discussed in this paper is applicable only to the 
M = 1 pure Yang-Mills, we strongly believe that there exists a renormalization scheme 
in which the supercurrent equations of Refs. [7J [11] which include matter fields are really 
justified. This "desirable renormalization scheme" should also justify various exact results 
in supersymmetric gauge theories. One of the properties which the desirable scheme 
should possess is the quantum action principle ( fl9|) . at least if we use the couplings which 
can be extended to background superfields. 

In the rest of this section, we would like to comment on how our result is consistent with 
the results obtained in Refs. [5j El [7J |8] . Grisaru, Milewski and Zanon [51 [6] studied the 
supercurrent equation by employing DRED. They constructed two supercurrents which we 
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will call J^di. an d i?ad i an d performed calculations at the two-loop level. One supercurrent 
J^2 satisfies the equation (in our notation and convention) 

D«j£> « D a (^^^[W a W a hn) , (42) 

where the coupling g and the operator [W^WqJdr are renormalized by minimal subtraction 
in DRED (i.e., the DR or DR scheme), and /3two-ioop is the two-loop beta function which 

(2) 

coincides with the NSVZ beta function at the two-loop level. The other supercurrent 
satisfies 

D<*J® « D a l^-^-^[W a W a }^ + bD 2 D a [vh^ (43) 

where /3 ne-ioo P is the one-loop beta function, and b is a coefficient which they did not 
explicitly compute. They claimed that jj^ contains the energy-momentum tensor, while 
the lowest component of satisfies the Adler-Bardeen theorem. 

Let us first consider Eq. f|42|) . At first sight, this equation may look inconsistent 
with our result, because it seems to give the trace anomaly which is not one-loop exact. 
However, this is not necessarily the case. The problem is that the normalization of the 
operator W a W a is quite ambiguous in DRED. We have discussed in section @] that this 
operator has an ambiguity of adding a term D 2 J, where J is a gauge invariant real vector 
superfield. However, in DRED, there exists the operator r^Ts (see section |3]) satisfying 

D 2 (m h ) = AeW a W a . (44) 

By choosing J oc (PT^), the normalization of the operator W a W a becomes completely 
ambiguous. This ambiguity may be related to the fact that there is no instanton in 4 — 2e 
dimensions, and hence the topological arguments do not work in this scheme. 

More concretely, suppose that the divergences of W a W a are of the form (1/ e n )Z) 2 (r' i r^) 
(n = 1, 2, • • •) [27], and we define a renormalized operator [W a W a ] as 

00 r (n) 1 

\w a w a ] = w a w a + J2^-jD 2 (m0 

71=1 6 4 

00 r (a) 

= (l + Cl (g))W a W a +^^-W a W a . (45) 

n=2 6 
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The n — 1 term eventually gives a finite counterterm. Then, we have a choice for the 
renormalization prescription. One choice is to set C\{g) = 0, which corresponds to the DR 
scheme. But we can also choose a scheme in which C\{g) ^ 0. The freedom in choosing 
c\(g) leads to the ambiguity of the renormalized operator [W a W a ]. 

One may consider that it is meaningless to say that the divergences are of the form 
(l/e n )D 2 (T^T0 instead of (l/e n )W a W a , since we have the relation flH}. In fact, this 
statement is quite nontrivial. The divergent pole (l/e n )D 2 (r /i r^) arises from n or higher 
loop diagrams, and hence the coefficient c n (g) is of the form c n (g) ~ g 2n + • • •. This 
means that the divergent pole (l/e n )W a W a comes from n + 1 or higher loops and has a 
coefficient of order g 2 ( n+1 ) . This is suppressed by additional g 2 than naively expected. In 
fact, this suppression can be confirmed by the two-loop calculation of Grisaru, Milewski 
and Zanon [5J, [6] . In our notation and convention, their result is summarized aslU. 

[W a W a ] DR = ( 1 - - e (tgpj + 0(g e )\ W a W a . (46) 

Having confidence that the divergence is indeed (1 / ' e n )D 2 {T^T , we suspect that there 
exists a "hidden divergence" of the form (l/e)D 2 (T^T0, which should be subtracted by 
nonzero C\{g). (In fact, Eq. f !2B|) is a similar kind of such hidden divergence which becomes 
clear by the analytic continuation of parameters into superspace [26i].) This subtraction 
will give a different normalization of the renormalized operator [W Q W a ] than the operator 
[r^ a ] DR subtracted by DR. 

The M = 4 regularization avoids the above difficulty by effectively replacing the 
counterterms as 

t C -^\D\Yh0 -+ D 2 K, (47) 

71=1 6 4 



where K is defined in Eq. (1351) . In this way, the structure of the divergence becomes 
much more clear. We can uniquely specify the overall normalization of [W a W a ], which is 
preferred by the topological nature of F^ V F^ U . (Strictly speaking, we have regularized the 



11 One should note that what they called the bare operator is in fact (g/go)W a in our notation, 
where go is the bare gauge coupling in DRED. The operator W a in their notation coincides with our W a 
because of the relation Zy = {g/go) 2 , where Zy is defined in their paper. The supercurrent should have 
been defined as {l/gl)W a W& + • • • instead of W^W^ + ■■■. 
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M = 4 theory by DRED to make our arguments concrete. But we believe that this is not 
essential. Any other regularization is possible as long as the arguments of section H] make 
sense.) Our renormalization prescription for W a W a gives the one-loop exact coefficient 
in the trace anomaly. 

(2) 

As to the second supercurrent J^a '■> our g uess is that this supercurrent is simply related 
to the first one by 

jM = ^±2PjW + [evanescent ] j (48) 

Ptwo— loop 

where by the term [evanescent], we mean an operator which vanishes in the limit e — > 
when it is inserted into a renormalized correlation function. In dimensional regularization 
or reduction, it often happens that two different-looking operators are in fact the same 
up to evanescent terms (see e.g., Ref. [32]). In fact, in the construction of the two super- 
currents, operators such as [D a , Da][T^T~\ DR are necessary which seem to be evanescent 
operators. We do not perform a detailed study on this point, but the existence of the two 
different supercurrents cannot be confirmed unless the possibility ( 148]) is excluded. (How- 
ever, the idea of two supercurrents becomes important when matter fields are included in 
the theory pH H21 E3].) 

Next we discuss the result obtained by Arkani-Hamed and Murayama [8J. They com- 
puted the anomalies in dilatation transformations. Their result can be summarized as 
follows. If we take the dilatation transformation of the gauge superfield V(x, 9, 6) as 

V(x, 6, 0) -> Vie'x, e t/2 6, e t/2 8), (49) 

where t is the parameter of the transformation, we obtain the anomaly given by 

J d 4 x J d 2 6 (t^^j W a W a + h.c. (50) 

On the other hand, we can canonically normalize the gauge field as V = g c V c , where g c is 
the "canonical gauge coupling" , and perform the dilatation transformation as 

g c V c (x, 9, 9) -> g'MVc^x, e t/2 6, e t/2 6), (51) 

where g' c {t) is chosen so that the gauge field is canonically normalized after the dilatation 
transformation. Then the anomaly is given by 

/ d*x J d 2 e\ - ^ log W*(g>m)W a (g>W c ) + h.c. (52) 
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Here g' c (t) is determined by the equation 



9?® 

Taking t infinitesimal, we get 



1 'it Ad 
L f _ 

9l 8^ 



tAd , 

871^ 



, 9c , 



(53) 



1 



9? it) 



J- , /3nsvz 

9l 9t 



(54) 



where /?nsvz is the NSVZ beta function. In the transformation ( 151 p . the classical action 
also changes. The total change of the path integral J[DV]e lS in this case is given by 





d 4 x J d 2 9t/3 NS y Z 



dgl 



Vc 



-^W a (g c V c )W a (g c V c )\ +h.c. 



(55) 



Therefore, it was claimed that there are two dilatation anomalies. One is given by Eq. ( 1501) . 
which has the one-loop exact form, while the other is given by Eq. ( J55l) which, at first 
slight, has higher order corrections. This result may seem to contradict with our result, 
since we have obtained the unique trace anomaly, up to a total derivative term and the 



ambiguity of choosing Ci discussed in Appendix A 



However, a closer look at Eq. (1551) gives us the resolution to this problem. We can 
rewrite this equation as 

/ OS 



p dS 

tpNSVZ 77-5- 

dg 2 c 



t/3 



NSVZ 



dgl 



+ 



2£6V , 



(56) 



where S is the action and 5S/5V is the functional differentiation of the action with respect 
to V. Naively, V(5S/5V) vanishes by equations of motion. However, the computation 
of the anomaly under the rescaling V — > (const.) V performed by Arkani-Hamed and 
Murayama indicates that the equation of motion in their regularization scheme should 
really be given as (see Appendix A especially Eqs. (1A.8[) . ( IA.9P and (lA.lOp ). 

1 



--(e 2V D a e~ 2V )W a ) + h.c. + (total derivative), (57) 



where one should note that the second term satisfies 

1 



1 

2 l 



:D *l-Z( e lV D a e -2V 



)W 



w a w a . 



(58) 
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Then, Eq. ( 1561) becomes 

t/W / dVO (-^ + J^L_) W«W a + h.C. + J d^xdWt^Ny 

= Jd^xJ d 2 6 (t^j W a W Q + h.c. + J d'xdVet^^Ny. (59) 

Therefore, the difference of the two anomalies is only the term containing Ny, which 
vanishes by the equation of motion. This just reflects the ambiguity of the definition 



of the energy-momentum tensor discussed in Appendix A In fact, by setting Sx^ = x M 
in Eqs. (IA.3[) and (I A . 5 P of Appendix A, we can see that the two dilatation anomalies 
computed above are just the trace T*} with different values of C/. In this way, we obtain 
the one-loop exact coefficient for the operator W a W a in the dilatation anomaly, which is 
consistent with our result. 

Shifman and Vainshtein [7] gave a supercurrent equation which is consistent with our 
result on the trace anomaly. Although there is no contradiction in the supercurrent equa- 
tion, their interpretation of the anomaly puzzle is different from ours. They argued that 
the one-loop contribution to the beta function comes from UV region, but higher loop 
effects are from IR region. Since operator equations are given in terms of UV quanti- 
ties, they claimed that the one-loop beta function should appear in the trace anomaly. 
However, as emphasized by Arkani-Hamed and Murayama [8], both the one-loop beta 
function and the NSVZ beta function can be derived by only using the information of 
UV physics. Indeed, the derivation of the NSVZ beta function we reviewed in section 2] 
has used only the decoupling of the adjoint chiral fields, Eq. ([28]) . which is completely 
determined in UV regiono Then the question of which beta function should appear in 
the trace anomaly comes back to us. This is precisely the problem we have studied in 
this paper. Although there are many renormalization schemes which give different beta 
functions, only one scheme possesses the desired property 

_d 

dg- 

for the appropriately renormalized operator [W Q W a ]#. Without this property, the trace 
anomaly formula is simply false. 



— 2 ((• ••)> = <((/ d*xd?6\[W«W a } R + h.c.) (• • •)) (60) 



2 We have not even used the Jacobian from the path- integral measure as in Ref. [5] . 
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Appendix A Energy-momentum tensor 

In this appendix we review basic properties of the energy-momentum tensor. The energy- 
momentum tensor can be defined as the Noether current of the Poincare symmetry. Let 
us consider a theory described by fields <fi{x) which are in some representations of the 
Lorentz group. We treat ipi as if they are bosons for simplicity, but a generalization to 
fermions is obvious. An infinitesimal transformation of the coordinates is given as 

x»^x'» = x> 1 + 5x l "(x). (A.l) 

We define a transformation of the field <pi under the above coordinate transformation as 

(pi(x) ->■ (p'j(x) = (fi(x) + 6ipi(x), (A.2) 
5<pi{x) = Sx^^ + ^Sx^cp^+d^dcpjix), (A.3) 

where Tlf is the Lorentz symmetry generator acting on </?/, and Ci is an arbitrary con- 
stant. This definition is taken so that it coincides with Poincare transformation when 5x^ 
is given as 5x^ = u^ u x u + a M , where and = —u VfM are constants. 
Let us consider a correlation function 

(n v^)) = / w^n^). (a.4) 

Then, by performing the above transformation of the fields in the path integral, we obtain 

°=-(V ^T^^n^(^o)+E(^/ ! (^)n^(^)) ) (a.5) 

\ i I i \ j=ii I 

where d is the space-time dimension. This is the defining equation of the energy-momentum 
tensor, but a slight ambiguity remains. There is a freedom to change T^ v as T^ u + d p A plll/ , 
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where A ppv is an arbitrary tensor satisfying A PPjV = —A ppv . We can use this freedom to 
make T^ v symmetric, T^ v = T UfJ- . First, note that if 5x^ = u p u x l/ where u pi/ is an antisym- 
metric constant tensor, the path integral measure times the factor e i5 ^' is invariant by 
the Lorentz symmetry Then / d d xT^LU pu = for an arbitrary constant uj pV) and hence 
is a total derivative, = d p B ppu . We can define a symmetric energy-momentum 
tensor as T pu -d p (B p ^ -B^ pu -B up ^), i.e., A ppv = -B p ^ v + B ppu + B vpp . In the following, 
we use this symmetric tensor. Even after making T pu symmetric, there still remains a 
freedom to change T pu as 

T pv -»■ T pv + d p d a C ppua , (A.6) 

where C ppua is an arbitrary field satisfying C ppua = —C piwa = C vapp . This is the so-called 
improvement of the energy-momentum tensor. If one prefers to define T pv by functional 
differentiation of the action with respect to the metric tensor g pu , this freedom corresponds 
to adding a term — t / (^R^py^C^ to the action, where R ppva is the Riemann tensor. 
Functionally differentiating Eq. ( 1A.5I) by 5x u (y), we get the Ward-Takahashi identity 

*E ((^ - vWvifa) ~ (\K + C h rT) (d y )^ d \ Xl - y) Vh {xS) JIM 3 *)) ■ 

(A.7) 

In particular, if y ^ Xi, this equation is just the conservation of the energy- momentum 
tensor, dpT^ w 0, where ~ means that the equation is valid up to contact terms. 

Aside from the improvement, there exists the ambiguity of choosing Cj in the definition 
of the energy- momentum tensor. To investigate this point, let us consider a rescaling 
transformation of the field (pj as 

(pj(x) (1 + a(z))^j(x), (A.8) 

where at(x) is an infinitesimal function. By performing this rescaling in the correlation 
function (1A.4|) and functionally differentiating the result by a(y), we obtain 



= (iN^Hip^Xi) ) + £ (* w 0* - v)v>/<(*i) Il^fe) ) • (A.9) 
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Here we have defined the operator iVj as 

Nj(x) = ^Uj(x) + Mx), (A.10) 

where Ai{x) is the anomaly arising from the path integral measure in the transforma- 
tion (IA.80 . The operator Nj vanishes by equations of motion, Nj ?a 0, since there are 
only contact terms in Eq. flA.9[) aside from the term containing Nj. Using Nj, one can 
see that the ambiguity in the definition of the energy-momentum tensor represented by 
Ci is of the form 

T» v = T fiV \ Cl =o-V^J2 C i N i- ( A - n ) 

i 

Thus, energy-momentum tensors with different values of Cj are the same up to the 
equations of motion, Nj ^ 0. More generally, any symmetric tensor which vanishes 
by equations of motion can be added to T^ v without violating the conservation equa- 
tion d^T^ fa 0. One can also check that this ambiguity does not affect the genera- 
tors of the Poincare symmetry P» = J d d ~ l xT^ and = J d d ~ 1 x(x l "T^ - x u T 0>M ) 
when Eq. (IA.7j) is integrated to give the commutation relations [P M , </?/] = id^ipi and 
[J^u, <pi] = i{xpd v - Xyd^ + E IfiV )(pj. 

Finally, let us comment on the role of Cj at the fixed point of RG flow. In section [2j 
we derive the formula for the trace of the energy- momentum tensor given by Eq. (j21|) . By 
using that formula, the dilatation current 

? D = x v TV-$ (A. 12) 

can be shown to satisfy the Ward-Takahashi identity at the fixed point d a \ a — (3 a = 0, 



£ ( t (d) & - y) {<d uVli { Xi ) + A It <p It (x t )) n <Pl 



jK Xj) 



+i(d y ),Y:(s id \x i - y) (-\y v T% ~ y'C^ tp^Xi) U rLi-r^ . (A.13) 

If we integrate this equation over to obtain the integrated version of the Ward- 
Takahashi identity, the final line of Eq. flA. 13[) does not contribute at all. Thus the 
arbitrary parameter Cj has no physical meaning in the scaling symmetry. 
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